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Abstract. We clarify and extend the theorem of Sveshnikov and Tkachov Jl[, ||, which gives an 
explicit connection between jet observables and the energy-momentum tensor. We check the relation 
between jet observables and the energy-momentum tensor for non-scalar (spinor and vector) fields, give 
a correct treatment of the light-cone singularity for massless particles, and extend the theorem of Q 
to the massless case. We also discuss the issue of gauge invariance. 



1 Introduction 

Modern QCD increasingly emphasizes precision measurements and perturbative calculations of higher 
order corrections (cf . measurements of a s and other parameters of the Standard Model) . 

It was argued in H that in the context of precision measurements a central role is played by a special 
class of observables — the so-called C-correlators — that contain all information about multijet structure 
and possess optimal stability properties with respect to experimental errors. As was shown in H, ||, the 
C-correlators possess another property that makes them extremely attractive from theoretical point of 
view. Namely, the theorem of Sveshnikov and Tkachov Q, Q expresses C-correlators (and, consequently, 
a vast class of other jet observables || ) in terms of the energy-momentum tensor in such a manner that 
no information about hadron bound states is used. 

In the arguments of Q, Q there are some gaps. The purpose of this work is to clarify them: 

(i) The theorem was proved in |Q, [Q only for scalar fields. So one needs to check it for non-scalar 
spinor and vector fields. 

(ii) The theorem was accurately proved only for massive particles. The massless case exhibits some 
subtleties due to the light-cone singularity that have to be clarified. 

(iii) The issue of gauge invariance in the case of gluons has to be clarified. 



2 Setup 

The C-correlators have the following form: 

(Y j ...Yj E ^'-- E iN /JV (PiL.-.P** ) ) • (1) 

\ ii in I p 

Here ik are indexes that run over all particles produced in an event, N is the order of correlator (N = 
1,2,.. .). Note that the definition is entirely in terms of observable quantities — particles' energies (Ei) 
and angles (pj). 

One can rewrite it in the Fock space formalism as follows: 

J dn^ . . . J dn N (m|e (ni)...E (njv)|m) X / w (ni,...njv), (2) 

where e (n^) is an operator- valued distribution on the unit sphere: 

£ ( n )^E/^ IPl a + (p)^(p)<5(P,n). (3) 

t Deceased. 



In fm, [2j the following relation between e (n) and the energy-momentum tensor was found: 
Theorem (Svcshnikov and Tkachov, 1995): 

1 

e (n) = lim t 3 / p 2 dp n t T lQ (pnt,t) , (4) 

t^oo J 


where Tjo is the energy-momentum tensor (the weak limit is implied here). 

In the context of high energy QCD we can regard gluon, quark and ghost fields as massless. This 
corresponds to high energy limit (p, n) = |p| ~ p . 

In order to compare the massive and massless cases, we will reverse the procedure of 0, 0] and from 
the field representation of T i0 obtain operator representation for both cases. 



3 Massive case 



3.1 Scalar field 



We already mentioned the result for scalar massive fields. Let us start from the well-known massive scalar 
Lagrangian and represent the fields in operator form: 



<P 0) 



(2tt 



Following [jij, Q one gets the following expression (in non-commutative case): 

oo 

/p 3 dp 
— — (a + a~(pn) + a~a + (pn)) . (6) 
ip 



3.2 Spinor field 

Now we turn to massive spinor fields for which we will obtain a similar result. We start with free 
asymptotic fields. Write 

, (l) = _L_/|L (e -,.„, (p)Wp)+e . f . l , (p)4(p) , m 

The expression for the energy- momentum tensor for free fields is well-known (see e.g. 0). So let us 
substitute (m) into (^) and use the stationary phase method. Only case of opposite signs in stationary 
phase equations contributes to e (n), because the same-sign contribution gives asymptoticaly e (n) = 0. 
By straightforward substitution (^) one obtains the following formula: 

l l 

lim t 3 [ p 2 dpn z ^- ^idoip - 9o^7i^] = lim — ^ / p 2 dpn l ( ( Tp- 
t^oo J 2 t^oo (2ttY J J 2p J 2q 



x exp{± " v ; " [(p) 2 - p 2 (n,p) 2 } T ° v ; P " [(q) 2 - p 2 (n,q) 2 ]} • F, (8) 
2m 2m 



F = |[U s (p')6+(p')7^ s (q')^(q') -^(p')d.(P , )7<«-(q / )^(q / )] 
Po 



where 



2 L ^(p>+( P ')7,^(q')^(q') -^(P')d.(p')7i^(q')4(q')]. (9) 

/ p , mp q 

r n +^i' q = /- (10) 
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The difference from the scalar case is only in the operator part of the expression. After rewriting e (n) 
in terms of 

P=|p| = ^= (11) 

and using for the operator part of expression (^) the following formulae: 

u s "ffj.u s = 2pfj,, Vs^/^Vs = 2p tJl , [d s ,df] + = (12) 



one obtains the final result: 



which agrees with the scalar case. 



o 



?^(bfb s (pn) + dfd s (pn)), (13) 



3.3 Vector field 

The massive vector field is treated similarly. One starts with the standard representation: 

M*) = / |V^(P) + H - C -)- ( 14 ) 

The vector field Lagrangian in general covariant arbitrary gauge is as follows: 

i 777 ^ CY 

L = -ts^A v a^A v + —A^A^ + -d»A»d„A^. (15) 

As a consequence of the Noether theorem it is always possible to add a 4-divergence to , which allows 
us to make a convenient choice of the Lagrangian. We starts with the Stuckclberg Lagrangian (a = 1 in 
(|l5|)) and use the fact that the theory in this case is gauge invariant. The Lorentz condition is applied 
(for a = 1 it appears automatically from the Lagrangian) in the operator form (which arises from fields 
structure). In this way one can avoid non-physical states. In this way there is no contradiction: one retains 
Lorentz invariance and local commutativity. But the components ^4^(p) are no longer independent. 

The negative-sign contributions to the Too can be eliminated by representing operators A^ (p) in ( fTi] ) 
in a local frame of reference in momentum space (cf. e.g. 0)). To obtain one can choose any gauge 
for example a = for calculational simplicity. The result is gauge independent and has the form of (|J), 
as expected. In the final formula only physical components survive: 

oo 

^( tar(pn) + a r a +(pn)). (16) 



4 Massless case 

The result ([!]) has to be interpreted carefully in the massless case. From the accurate derivation (given 
below) it follows that the integration from to 1 should be spread (formally integration with 6*-function 
over p) over region A (see Fig.l on the next page). 

At infinitesimal p we can neglect contribution of particles in the A— region, because of limited 
experimental sensitivity and hence such slow particles could not be registered by detectors. In the 
massive case, the region Ai simply does not contribute as follows from the formula for stationary phase 
manifold: in the denominator square root of negative value does not have any physical meaning. But 
presence of the region Ai plays an important role when we discuss massless fields. 



3 



*(p)a 



fl(i-p) A 



1 p 

Fig.l 



1 p o 




i p 



Fig.2 



Fig.3 



4.1 Scalar and spinor fields 

Consider the massless case. Here situation becomes rather intricate. Integration over interval [0,1] in ([|) 
- if done formally — leads to a result that differs by coefficient 1/2 from what is to be expected from 
analogy with the massive case. The difficulty is due to the fact that all massless states are sitting on 
the light cone so that their distribution $(p) = p 2 riiTio(pnt,t) is a 6— function. However, the position 
of the latter corresponds to the boundary of the integration region at p = 1. Therefore, care is needed 
in how one defines the S— function. The trick we use ensures a natural connection of the massive and 
massless cases. Namely, let us extend the integration region into positive direction over unit interval to 
1 + e and decrease it near zero, obtain: p = [1 — e' , 1 + e] ,Ve' <E (0, 1) , Ve € (0, oo). The non-zero part 
of the integral ([|) comes from the infinitesimal region near 1 symmetrically on both sides. 

All massless particles occupy the infinitesimally narrow <5-region at p = 1. The stationary phase 
method gives the equations: 

\V\ =pn ' W\ =pn ' (7) 

e ■ e'np + nq = 0, (18) 

where e and e' independently takes the values ±1. For the positive sign of e • e' , the stationary manifold 
is p = q = 0, |p| = po = 0. In the stationary region asymptotically e (n) = 0. 

For the negative e • e' sign the stationary manifold is: p = 1, p = q = tun, in the stationary region 

p' = zun+-^, q' = tun+-^=, p ' = \-JL. (19) 
yt \Jt \Jt 

Here after representing p, q as a series we use change of variables p <-> p',q^> q'. As was announced 
above, let us introduce e into integration limits over p (0 < e < 1). After substituting the above into (0), 
for scalar fields one obtains the following formula: 

e (n) = km ^ / p' Hp' * f # J g,' ■ ^ • e • e> • .'^W) 

e 

x exp {it [e (p' - p'np') + e' (<Zo - q'np')]} (20) 

oo 

Let us insert into this expression 1 = J S(tu — p) dm. One can represent the absolute value of p as series 

o 

and retain only contributions of order 0(^=). It yields the argument of the J-function: ^— j ■ One 
must take into account 

(p x ,n) = 0, p = p±+p,|. (21) 

The integrals over p_i_,q_i_ are trivial and yield the factor (2tt) 2 w 2 . Now one can do the integrals over 
qu , p in two independent ways with the same final result: 
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The first way is to take the integral over p first and apply the formula: 



dqn^. >- = 2-k. (22) 

% % 

— oo 

Note that the result is independent of a, b > 0. 

The second way is to take the integral over q» first and then the integral over p. Take the limit and 
use the formula: 



dpS(p) = 1. (23) 
In both cases one obtains the following result: 

oo 

Ifrffola+a-W + a-a+trn)). (24) 
o 

For massless spinors the calculations are similar. The result is 

oo 

- J m 2 dw (b+b s (mn) + dfd s (wn)) . (25) 



4.2 Gauge field 

The massless vector case is slightly more interesting because it is necessary to consider the problem due 
to broken Lorentz invariance when one deals with physical states. It is impossible to use the local frame 
of reference in momentum space as in the massive case because of the denominator that is singular at ( 
m = ) . To circumvent the problem one can use the standard Gupta-Bleuler quantization procedure (see 
e.g. Q). For quantum fields the Lorentz condition can be rewritten as a condition for physical states: 

d ll A-$ = Q, <Z>*d„A+ = 0, (26) 

Now one turns to the local frame of reference || . The non-physical zeroth and third components cancel 
each other as expected. In the average value over a physical state $ only the physical components 
remain. (Average value of observable over physical states is the same as over all states). Let us turn 
to the Lagrangian in Lorentz gauge (a — 0). One takes into account the above arguments, repeats the 
calculation of e (n), and obtains: 

oo 

o 



5 Conclusion 

We have seen that the expression ([!]) of C-correlators in terms of the energy-momenum tensor is 
correct for non-scalar particles as well. The distribution of massive particles lies inside the integration 
region for p (see Fig. I on the previous page). 

In the massless case the essence of calculation is that all particles are sitting on the light cone. 
Therefore the operator- valued distribution $ is <5-function (see Fig. 2 above). It follows from our analysis 
that the ^-function that describes the integration region must be defined as a limit "from the right" of 
smooth regulators (Fig. 3). 

This definition for the massless case should be important: (i) in axiomatic proofs; (ii) in proofs of a 
generalization of KLN theorem for jet observables Q 

*This work was planned by N.A. Sveshnikov. 
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For gauge particles we verified the formula (^) checked for Lorentz gauge using in the Gupta-Bleulcr 
formalism. It would be interesting to develop a treatment for general (not necessarily covairant) gauges. 
There is systematic formalism for treatment of arbitrary gauges for the case of Green's functions (see 
e.g. ||) but not for the Fock space. It would be useful to develop such a formalism in order to clarify 
calculational issues enocountered when there are loop and phase space integrals simultaneously. 
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